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We study the confidence and backaction of state reconstruction based on a continuous weak 
measurement and the quantum filter equation. As a physical example we use the traditional model 
of a double quantum dot being continuously monitored by a quantum point contact. We examine 
the confidence of the estimate of a state constructed from the measurement record, and the effect 
of backaction of that measurement on that state. Finally, in the case of general measurements we 
show that using the relative entropy as a measure of confidence allows us to define the lower bound 
on the confidence as a type of quantum discord. 

PACS numbers: 03.65.Ta, 03.65.Yz, 73.63.Kv 



I. INTRODUCTION 

Recently there has been a great deal of activity on 
the topic of "weak" quantum measurements [lj-|7 | in both 
mesoscopic (8l-[Toj and macroscopic systems [11 4l3| . In 
contrast to projective measurement, weak measurements 
only perturb the system slightly, but in turn can only 
provide limited information. According to the theory of 
open quantum systems, both the evolution of the quan- 
tum state and its decoherence depend on the system- 
apparatus coupling strength and the basis in which the 
measurement system is measured. One of the advantages 
of such a measurement is that, given a weak continuous 
measurement record one can reconstruct the quantum 
system state during its evolution. One particularly inter- 
esting approach, which we apply and investigate here, is 
the "quantum filter equation" as pioneered in the quan- 
tum limit by Belavkin [l5| and others [l6j . 

The quantum filter equation has been shown to be a 
powerful method for state reconstruction, and is fairly ro- 
bust in terms of the resolution needed in describing the 
measurement record. For example, in recent work [17( it 
has been shown that the continuous "analogue" measure- 
ment record can be reduced to a "one-bit record" and still 
the filter equation can efficiently produce an estimate (or 
purification) of the system state. Similarly it has been 
shown that feedback control can be used to enhance the 
speed of the state estimation [l8l - l20j , and that it can be 
further optimized when combined with a kind of process 
tomography [2l|. However, successful application of the 
filter equation requires accurate knowledge of the evolu- 
tion, both coherent and incoherent, that the monitored 
system undergoes [22j . 

Here we investigate how this state estimation method 
can be optimized by considering both how confidently an 
estimated state deduced from the measurement record 
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reflects the actual state of the system, and how the effect 
of backaction changes the state. Here the "actual state of 
the system" means the system under the influence of the 
back-action [23-25] of the measurement apparatus, not 
the initial prepared state. In other words, we quantify 
the confidence of the state-estimation process indepen- 
dently from the overall fidelity of the measurement. With 
these two quantities, confidence and backaction, one can 
identify the most sensitive parameter regime for estimat- 
ing the original system state. Such quantities have been 
employed in investigating information gain with projec- 
tive and general measurements [H, [53] • The behavior of 
these quantities in the context of weak continuous mea- 
surement has not been studied in great detail as of yet, 
though the concept of informati on g ain and measurement 
disturbance is well understood [Hf [2t| . In addition, we 
point out that our approach here is different from that 
used in some other works. For example, here we are con- 
cerned primarily with the optimal measurement of an un- 
known state by refining our state of knowledge, whereas 
in some other approaches the goal is primarily manipu- 
lating (or purifying) one's state of knowledge of a given 
quantum system, and the initial unknown quantum state 
is unimportant [29(. 

The model we use here is based on the continuous mea- 
surement of a double quantum dot (POP) charge state 
using a quantum point contact (QPC) [8|, [3fj-[32j . In this 
situation the tunnelling barrier of the QPC is modulated 
by the charge in the nearby DQD, and produces a con- 
tinuous measurement record which can be used in the 
filter equation. However, the approach is quite general, 
and recent applications have also arisen in circuit-QED 
systems [UHj]- This article is organized as follows: In 
Sec. II we provide a general definition of the confidence 
and backaction. In Sec. Ill we introduce the model for 
weak measurement of a DQD, and show the results given 
by the quantum filter equation. In Sec. IV we specu- 
late about a possible filter equation based on ensemble- 
averaged measurements. Finally, in Sec. V we show that 
using the relative entropy as a measure of confidence al- 
lows us to define the lower bound on the confidence as 
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FIG. 1: (Color online) (a) Schematic of the components of the 
quantum state estimation using the filter equation. Here, p s 
is the system initial state, pjj(t) is the evolution of the initial 
state of the system in the measurement-induced noises, pfi(t) 
is the estimation of the system state, and T>[pE, pn\ is the esti- 
mation confidence, (b) shows a diagrammatic representation 
of the confidence C and backaction B as the distance between 
the various states. Note that when we use the relative en- 
tropy as a distance measure these quantities are asymmetric 
(as represented by the one-way arrows). 



the quantum discord. 



II. DEFINITION OF CONFIDENCE AND 
BACKACTION 

Consider a quantum system S which interacts with a 
measurement apparatus A. When implementing a gen- 
eral quantum measurement M. = {H}, our knowledge of 
the system state is based only on the measured appara- 
tus output y(t). We denote our estimate of the system 
state as pE{t). The initial state of the system is denned 
as p s (0), and p R {t) is defined as the evolution of the 
initial state of the system in the measurement induced 
noise, i.e., what we can think of as the actual state of 
the system. We define the ideal state of the system, i.e., 
its evolution if it were not affected by the measurement 
apparatus at all, as pi{t), i.e., the entirely coherent evo- 
lution of p s (0)- 



Definition: The quantum state estimation confidence 
is defined as 



C = D[p E ,pn], 



(1) 



where D[-\ is some appropriate distance measure. Gener- 
ally speaking, the smaller C is the more confident we are 
of the estimated state, and C = if and only if pe = Pr, 
which means that the estimated state is totally confident. 
Similarly we define the backaction by 



B = D[p!,p R ], 



(2) 



so that B — implies no backaction. In the treatment 
that follows of course we have full knowledge of all these 
states at all times, and can thus in a theoretical sense 
identify the optimal parameters that minimize both of 
these quantities. 

There is some freedom in choosing an appropriate mea- 
sure for both C and B. Here we explicitly consider both 
the fidelity and the relative entropy [34[. The fidelity is 
commonly used to measure the effectiveness of the filter 
equation, and is defined as 



F = 1—C 



(3) 



We define the confidence as C = 1 — F in this case, to 
match our definition of a distance measure, so that high 
fidelity implies C = 0. However since the fidelity is a 
pseudo-distance this lacks some characteristics of a true 
measure. In the case of the relative entropy we define, 



where 



C = S(p R \\p E ), B = S( Pl \\p R ), 



S(a\\p) = -Tra log p- S(a) 



(4) 



(5) 



The relative entropy can be seen as a distance measure, 
though as it is asymmetric in a O p it is technically not. 
In fact the ordering here is important; with the inverse 
ordering the backaction B — » oo as pi is a pure state 
in this example. Using the relative entropy allows us 
to make a more direct connection to a general positive 
operator valued measure (POVM) description of a weak 
continuous measurement in the final section of this work. 



III. CONTINUOUS WEAK MEASUREMENT 
OF A DOUBLE QUANTUM DOT 

We now present the specific details of the DQD and 
QPC system. A DQD consists of a dot L, connected to an 
emitter, and dot R, connected to a collector |30l. l35l. l36|. 
As is typical, we assume that the DQD is in the strong 
Coulomb regime such that only one electron is allowed in 
the whole DQD. Here we assume the DQD is prepared 
in a single electron state, then isolated from the emitter 
and collector reservoirs via manipulation of gate voltages. 
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The two single-dot states are denoted by \L), and \R). 
The Hamiltonian of the DQD can be written by 



Here 



Hdqd = 0<7 x /2 + ea z 
\L)(R\ + \R){L\, 



(6) 
\R)(R\, 



with a x 

ujq = \J t 2 + 4|£1| 2 , e is the level splitting between the 
two single-dot states and Q is the coherent tunnelling 
amplitude betweens the two dots. 

In this simple model, the current shot noise |37| of the 
QPC acts as a noise source. In the low-temperature limit 
with ksT <C foui, the noise spectral density takes the form 

MM 



Rk 



D(l - D) 



(eVqpc - fiw) 



[1 - exp{-(eVQ PC - huj)/k B T} 



,(7) 



where D is the transparency of a single channel in the 
QPC, and Rk = h/e 2 is the von Klitzing constant. To 
treat the measurement signal, or current through the 
QPC, as a classical variable one must assume that the 
QPC evolution is much faster than the DQD, so that only 
the zero-frequency component is important; this is effec- 
tively a Markovian limit in terms of treating the QPC 
backaction on the DQD [jl^. 

We also treat the QPC as a perfect detector. In this 
limit we can define the measurement strength of the QPC 



(A/) 2 
16J(0)' 



(8) 



Here, AI = I L — I R , where I L = D L e 2 V / irh,I r — 
D R e 2 V/irh, D L = D + AD, D R = D - AD, and AD 
is the change in the transmission of the QPC due to the 
charge state of the DQD. 

The evolution of the real state of the DQD p R can be 
derived using the Bayesian techniques of Korotkov [8Hl"oj 
to give the selective stochastic master equation (SME) in 
Ito form, 



dpR 



[Hbqu-.Pr] dt + kT>{<j z } p R dt 



+ V2^n[a z ] PR dW R , 
where k is defined above, 

rt 



(9) 



H[a z ]pR = a z p R + p R o\ - Tr(a z p R + p R al)p R , (10) 

and the real Wiener process dWR satisfies E(cIWr) = 0, 
(dWR) 2 = dt. In Eq. the super-operator V is defined 
as 



T>[a]p = apa) — — (a)ap + pa) a). 



(11) 



Given that, in a general sense, the measurement oper- 
ator is y — V 2nha z the measurement record increment 
at a time t is, 



dy(t) 



^{a*{t))dt + dW R = AJ j <7z dt + dW R . 

V 2J(0) 

(12) 



(<r? (t)> = Tr[<T zPR (t)} 



(13) 



is the instantaneous expectation value of a z at time t 
based on the selectively evolved density matrix Pr{€). 

This equation of motion also relies on the assumption 
that 



\I L - Ir\ = \AI\ « I = (I L + I R )/2, 



(14) 



so that many electrons, N > (Iq/ AI) 2 3> 1, should pass 
through the QPC before one can distinguish the quantum 
state. This is the weakly responding or linear regime, and 
the model as we have described it is entirely equivalent 
to the formulation used by Korotkov and others [8T4l0l|. 
Also, note that for consistency with other works on the 
filter equation as a state estimation technique [13, HH we 
describe the noise as a Wienerprocess, so that the width 
of the Gaussian distribution 0, [Io| used to describe the 
weak measurement with a QPC is absorbed into K. 



A. Quantum Filter Equation 

To estimate the quantum state pe from the mea- 
surement output we employ the quantum filter equation 
method {[Mm IS E|- The evolution of the estimated 
state P E is described by the following stochastic mas- 
ter equation, identical to the "system" one, except the 
measurement signal from the system evolution, described 
above, determines the noise term: 



dp B 



[-HdqDi Pe] dt + KV[a z ]p E dt 



1kH.\o z \pe 



AI 



y/2J{0) 



{(^(t))-(af(t))}dt + dW R 



(15) 



The last term is analogous to the classical innovation 
in control theory [10, i.e., the difference between the 
actually measured current and the predicted current with 
the estimated state. The state estimation process in- 
volves setting P e(0) — 1/2, and evolving under the noise 
determined by the measurement record from the "exper- 
iment", or in this theoretical work Eq. ©. 

As demonstrated in Ref. [22| . even a small error in 
the Hamiltonian of the above equation can induce errors 
in the estimate of the state provided by the quantum 
filter equation. We expect the same to be true of the 
estimates of the noise spectrum of the QPC. Here our 
goal is to study the efficiency and robustness of quantum 
state estimation via the filter equation, so as in Ralph 
et at (171 we choose the same Hamiltonian and QPC 
properties in both Eq. (fTB"]) and Eq. ©. We leave the 
problem of accurately estimating the Hamiltonian and 
the noise properties of the measurement in a dynamic 
way [22j for future work. 
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FIG. 2: (Color online). The top row of figures shows the ensemble-averaged behavior (over 2,000 realizations) of (a) the 
confidence (blue) and backaction (red) defined in terms of the fidelity, (b) the occupation of the left "dot" or state for the 
system, Eq. ([9]) (in red) and the estimator, Eq. fj 1 5j) . (in blue), and (c) the confidence (blue) and backaction (red) defined in 
terms of the relative entropy. In all figures k = 0.005O, and time is evolved for the equivalent of 15 Rabi oscillations of the bare 
quantum dot state. Figure (b) shows more directly how the estimator quickly oscillates in phase with the system state, but 
takes time to evolve to the same population magnitude. The small oscillations seen in both backaction curves in (a) and (c) 
are typical for the definition of the backaction, and simply represent the "nodes" in the oscillation curve of, e.g., Pl, where the 
pi and pr states coincide. Figures (d)-(e) show the same quantities as (a)-(c) except just a single realization. The estimator 
state still quickly approaches the system state, as is typical with the filter equation. The backaction changes in magnitude 
drastically in both (d) and (f). There is no single-realization dephasing in (e) because we assume the QPC to be a perfect 
detector. This is in contrast to recent work [13 ] on circuit QED where some information can remain hidden due to the finite 
lifetime of the measurement cavity. 



Finally, we combine the time-evolution of these two 
equations © and (|T5|) to calculate the confidence of 
the estimated state and backaction of the measurement 
using both the fidelity and the relative entropy, as de- 
scribed earlier. In Fig. 2, we show numerical results for 
these quantities. We set pr(0) = \L)(L\, /j.e(0) = 1/2, 
and evolve using the standard techniques for 150000 time 
steps per Rabi oscillation. The top row of figures shows 
the results averaged over 2000 realizations, while the bot- 
tom row shows just a single realization. In these results 
we set h = 1, e = , f2 = 1 and k = 0.005J7. Com- 
paring to real parameters from [3fjl |39| , one could choose 
a strong inter-dot coupling of the order of fi = 32 /zeV, 
giving a timescale of 130 ps for the Rabi oscillations we 



show in Fig. 2. This should be chosen carefully however, 
to match the desired properties of the QPC timescales 
(or whatever the measurement device happens to be). 

When we acquire information from the measurement, 
it of course induces significant backaction on the system 
itself. Figure 2. (a) shows the confidence and backaction 
in terms of the fidelity, while Fig. 2.(c) shows the same 
in terms of the relative entropy. Both give reasonable 
descriptions of the distance between the estimated state 
and system state, and for weak measurement strength 
the confidence saturates before the backaction does. Ob- 
viously then the trade off is to measure on a time scale 
where both the confidence is relatively high and the back- 
action is low. 
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To gain more insight on what is actually happening 
during the evolution of the filter equation, Fig. 2(b) and 
(e) show the population of the left state of the dot for 
both the system and estimator. In the ensemble aver- 
aged case Fig. 2(b) backaction from the measurement 
dephases the state, but the estimator matches the sys- 
tem state before coherent information is totally lost. In 
a single realization, 2(e), the system state does not de- 
phase because the QPC acts as a perfect detector. Com- 
pare this to the case of circuit-QED where the lifetime 
of the cavity can induce dephasing on certain timescales 
[lil HI)]. The off-diagonal elements behave in a similar 
fashion. 



justify this correspondence, though one can note that Eq. 
(fl6|) represents a valid solution to Eq. (fT5j) for the tra- 
jectory determined by dWn = 0. We also point out that 
the fictitious non-linear force in the second line of Eq. 
(|16p is not physical, and the equation may not ensure 
positivity of the density matrix pE at an arbitrary time. 
Why this works so well in reproducing results from the 
ensemble averaged filter equation, at least for this case 
of a single qubit measured in the a z basis, is not clear, 
and represents a possible avenue of future work. 

V. POVM AND THE DISCORD AS A BOUND 
ON CONFIDENCE 



IV. A FILTER EQUATION FOR ENSEMBLE 
EXPECTATION VALUES 

Solving for the ensemble-averaged results by collat- 
ing many single realizations is sometimes an arduous 
numerical task, though can be useful in the stochastic 
Schrodinger form if one is modelling a system with a large 
Hilbcrt space. In practise the system state density matrix 
ensemble averaged over all measurement trajectories can 
be trivially calculated by averaging Eq. ([9]), and noting 
.EfdWfl] = 0. This gives the expected Lindblad equation 
of motion which induces the behavior we observe in the 
backaction. How about the estimator state? Averaging 
Eq. (|15p is non-trivial as the individual trajectories de- 
termined by (crf'(t)) are not statistically independent of 
PE(t). 

Rather than attempt to do so we simply write down 
an analogy to the quantum filter equation which de- 
pends on ensemble-averaged quantities rather than in- 
dividual trajectories. We now define pe — E[pg 
{af{t)) = E[(af(t))], and (af (*)) - E[(af (*))]. Thus 
the term Fi[(af (t))] represents the expectation value ex- 
tracted from an ensemble averaged version of Eq. ©, i.e., 
the evolution of the real system under the effect of de- 
phasing. By comparison to the stochastic filter equation 
we consider the following non-stochastic filter equation, 



^ = -t [Hdqd,Pe] + nV[a z ]p E 



+ V2kH[<t z ]pe 



Al 



(16) 



Solving this equation directly is computationally trivial. 
To illustrate this we plot the confidence as a function of 
time and measurement strength n in Fig. 3. We can 
easily see that as k increases the confidence saturates 
quickly, but with an associated strong backaction, as ex- 
pected. Remarkably, if we inspect the density matrix 
elements of the estimated state generated by Eq. (fT6| to 
those generated by the ensemble average over trajectories 
of Eq. (|T5j) they coincide closely. This is illustrated in 
Fig. 3(c) and (d). Curiously we are unable to rigorously 



In quantum theory one can describe any measure- 
ment scenario as a positive-operator valued measure 
(POVM). For example, weak measurement is some- 
times discussed in terms of a POVM on a combined 
system/measurement-apparatus, where the measurement 
apparatus itself is also considered to be a quantum sys- 
tem. To gain a more fundamental perspective on the 
confidence and backaction, as we have defined them here, 
we consider an alternative measurement scheme of an 
asymptotic POVM. 

First, we retain our definition of the initial pure 
system-apparatus state pi(0) = p s (0) ® Pa(0). We then 
assume that a measurement apparatus A is allowed to in- 
teract with the system to produce the typically entangled 
and correlated system-apparatus state p Si A (we suppress 
any time argument here for complete generality). We de- 
fine the "real" state of the system then as this combined 
state p s ,A- Finally, we perform measurements on A asso- 



ciated with a POVM {ir/' T IP/}, where £V Ilf' r Tlf = 1. 
J j Our estimate of the combined system-apparatus state 
given by the measurement is p m = II^p.bII- . Again 
we define the confidence in terms of the relative entropy, 
so C = S(p a<A \\p m ). 

The relative-entropy-based confidence has an interest- 
ing lower bound in the case when the POVM becomes a 
projective valued measure. Writing, 



c = -Tr Ps , A iog^n/p 5 , A nf 1 - s( Ps 



(17) 



we can substitute lif = where \j) is some distin- 

guishable orthonormal basis describing the measurement 
apparatus. Then the confidence becomes, 

C = -TrY^\j){j\p s ,AlogJ2\3){j\ps,A\j){j\-S(p s ,A) 
j j 

= -Tr J2 U> (J\Ps,A |i) (j I log ]T | j) (j\p s ,A\j) (J I 



-S(p s ^A) = S(pm) - S(p S} a) 

> min| j )[S'(p m )] - S(p a> A) = V, 



(18) 



where T> is the discord as defined by Modi et al [34| , when 
they assume classicality in terms of only one sub-system. 
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FIG. 3: (Color online) Ensemble-averaged (a) confidence C 
and (b) backaction B in terms of the relative entropy between 
system and estimator states, as a function of the interaction 
time t and measurement strength k. Both figures are derived 
using the ensemble-averaged equation of motion ()16[) . (c) and 
(d) show a comparison between the probability of occupation 
of the left dot Ph(t) derived from averaging Eq. [15] over 
many realizations (dashed blue lines) to that given by solving 

Eq. nnj. 



In the case of a general POVM, one could argue that the 
minimization of C over all possible POVMs is equivalent 
to a generalization of the definition of Modi's discord. 
Finally, we note that there is a correspondence between 
the estimator state pe we discussed in terms of the filter 
equation, and the partial trace Tr^(/o m ) over the state 
constructed from asymptotic POVM measurements (and 
the same for the real state pn and Tta(Ps,a) evolved in 
the measurement noise in the filter equation example). 



VI. CONCLUSION 



In many realistic quantum readout architectures the 
reliability of the quantum measurement output is an im- 
portant issue. In this article we discussed how to mea- 
sure the confidence and the backaction of a state recon- 
structed from continuous weak quantum measurement. 
As a typical example, we considered a DQD measured 
by a nearby QPC. Based on the theory of open quan- 
tum systems and the quantum filter equation method we 
briefly discussed the trade-off between measurement con- 
fidence and measurement-induced backaction. We also 
considered a possible filter equation for ensemble aver- 
aged results. We finished by discussing the case of a 
general POVM and how the confidence (when defined as 
a relative entropy) has a lower bound related to Modi's 
discord [341 ] . 
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